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ABSTRACT 


During  the  summer  of  1999,  the  Applied  Physics  Laboratory  worked  with  the 
Alaska  Department  of  Fish  and  Game,  Biosonics  Inc.,  and  Hydroacoustics  Technology 
Inc.  to  make  underwater  environmental  acoustic  measurements  in  Alaska’s  Kenai  and 
Wood  rivers.  The  goal  of  these  measurements  was  to  understand,  and  model,  the 
sources  and  magnitude  of  background  acoustic  reverberation,  defined  as  the  average 
acoustic  intensity  versus  time  (range)  in  the  absence  of  any  scattering  from  fish. 
The  background  reverberation  determines,  in  part,  the  minimum  size  class  of  fish 
detectable  and  the  accuracy  of  fish-tracking  and  counting  algorithms  that  exploit  the 
complex  phase  of  the  fish  echo.  It  is  also  necessary  to  establish  rigorous  detection 
probabilities  for  all  fish  size  classes.  This  report  summarizes  the  results  and  analysis 
of  these  measurements  and  presents  environmental  acoustic  models  that  can  be  used 
to  evaluate  signal-processing  and  fish-tracking  algorithms. 


TR  2001  iii 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


CONTENTS 

EXECUTIVE  SUMMARY  .  vii 

1.  INTRODUCTION  .  1 

2.  MEASUREMENTS  IN  THE  KENAI  RIVER  .  3 

2.1  Site  Description  and  Geometry  .  3 

2.2  Acoustic  Measurements  .  5 

2.3  Reverberation  Due  to  Volume  Scattering  .  7 

2.4  Reverberation  Due  to  Bottom  Backscattering  . . .  12 

3.  MEASUREMENTS  IN  THE  WOOD  RIVER  .  16 

3.1  Site  Description  and  Geometry  .  16 

3.2  Acoustic  Measurements  . 17 

4.  SUMMARY  AND  CONCLUSIONS  .  22 

REFERENCES  .  25 

APPENDIX,  Comparative  Measures  of  Scattering  from  Bubbles 

and  Glacial  Silt  Pau’ticles  emd  a  Model  Relating  Turbidity 

Measures  to  Volume  Scattering  from  Glacial  Silt  Particles  .  A1 


iv  TR2001 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


LIST  OF  FIGURES 


Figure  1.  Map  of  the  lower  Kenai  River  showing  location  of  the  sonar 

test  site  .  3 

Figure  2.  Experimental  geometry  for  acoustic  reverberation  measurements 

made  in  the  Kenai  River  and  key  geometric  variables  .  4 

Figure  3.  Results  of  the  CTD  measurements  made  in  the  Kenai  River 

at  1030  local  time  on  8  June  1999  .  5 

Figure  4.  Image  of  100  pings  transmitted  at  1-s  intervals  and  a  frequency 
of  420  kHz  versus  ping  number  and  range  from  sonar; 
reverberation  level  from  ping  number  33,  mean  RL  based  on 
the  100-ping  ensemble  average,  and  a  model  for  RL  .  9 

Figure  5.  Mean  RL  at  200  kHz  versus  range  from  sonar  based  on  a 

100-ping  ensemble  average  and  a  model  for  RL  .  10 

Figure  6.  Sample  histogram  of  the  ping-by-ping  420-kHz  reverberation  level 

corresponding  to  a  range  of  15  m  in  the  data  shown  in  Figure  4  . .  11 

Figure  7.  The  geometry  of  backscattering  from  a  sloping  riverbed  .  13 

Figure  8.  Mean  RL  measured  at  200  kHz  versus  range  from  sonar 
based  on  100-ping  ensemble  average;  model  for  RLy] 
model  for  and  the  incoherent  sum  of  RLy  and  RLi,  .  14 

Figure  9.  Mean  RL  measured  at  420  kHz  versus  range  from  sonar 
based  on  100-ping  ensemble  average;  model  for  RLy] 
model  for  RLi,,  and  the  incoherent  sum  of  RLy  and  RLt,  . .  15 

Figure  10.  Map  of  the  Wood  River  showing  the  location  of  the  sonar 

test  site  .  16 

Figure  11.  Experimental  geometry  for  the  acoustic  reverberation 

measurements  made  in  the  Wood  River  and  key  geometric 
variables  .  17 

Figure  12.  Mean  RL  measured  at  200  kHz  versus  range  based  on  a 

100-ping  ensemble  average,  and  a  model  for  RL  .  18 

Figure  13.  Mean  RL  measured  at  420  kHz  versus  range  based  on  a 
100-ping  ensemble  average  for  9p  —  0°  and  6p  =  8°, 
and  a  model  for  RL  .  19 


TR  2001  V 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


Figure  14.  Sample  histogram  of  the  ping-by-ping  420-kHz  reverberation 
level  corresponding  to  the  data  shown  at  a  range  of  5  m 
in  Figure  13  . .  20 

Figure  Al.  Backscattering  cross  section,  for  a  single  bubble  or  particle 
of  glacial  silt  versus  bubble  or  particle  radius,  a,  for  acoustic 
frequencies  of  420  kHz  and  200  kHz  .  Al 


vi  TR2001 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


EXECUTIVE  SUMMARY 


During  the  summer  of  1999,  the  Applied  Physics  Laboratory,  along  with  the  Alaska 
Department  of  Fish  and  Game,  Biosonics  Inc.,  and  Hydroacoustics  Technology  Inc., 
conducted  sonar  measurements  in  Alaska’s  Kenai  and  Wood  rivers.  The  goal  of  these 
measurements  was  to  gain  an  understanding  of,  and  model,  the  sources  and  magnitude 
of  background  acoustic  reverberation.  Background  reverberation  is  defined  as  the 
average  acoustic  intensity  versus  time  (range)  in  the  absence  of  any  scattering  from 
fish.  The  background  reverberation  determines,  in  part,  the  minimum  size  class  of 
fish  detectable  and  the  accuracy  of  fish-tracking  and  counting  algorithms  that  exploit 
the  complex  phase  of  the  fish  echo.  In  addition,  it  is  necessary  to  establish  rigorous 
detection  probabilities  for  fish  of  all  size  classes.  The  measurements  were  at  200  kHz 
and  420  kHz,  in  conditions  and  geometries  that  simulated  the  configuration  used 
by  sonars  for  counting  fish  in  rivers.  These  measurements  produced  estimates  of  the 
reverberation  level  (RL)  (in  dB  re  fjFa)  and  estimates  of  the  variance  of  this  quantity. 

There  are  two  primary  sources  of  reverberation,  scattering  by  particles  and  bubbles 
entrained  in  the  water  volume  {volume  scattering)  and  scattering  from  the  riverbed 
{bottom  scattering).  Our  study  shows  that  volume  scattering  is  the  more  important 
of  these  two  contributions,  and  it  can  be  described  by  the  scattering  strength  (5„)  in 
decibels. 

For  the  Kenai  River  site,  we  estimated  Sy  to  be  —67  dB  ±4  dB  at  200  kHz  and 
-53  dB  ±4  dB  at  420  kHz.  The  difference  is  predicted  by  the  frequency-to-the-fourth- 
power  law  associated  with  Rayleigh  scattering.  For  the  Wood  River  site,  Rayleigh 
scattering  also  applies,  and  we  estimate  Sy  to  range  from  -55  dB  to  —66  dB  at 
200  kHz  and  from  -42  dB  to  -53  dB  at  420  kHz;  this  larger  range  is  associated  with 
tidally  influenced  river  currents  which  modulate  the  bubble  density.  We  also  modeled 
fluctuations  in  RL  and  demonstrated  that  at  a  fixed  range  and  under  stationary 
conditions,  fluctuations  in  RL  and  Sy  are  well  described  by  the  probability  density 
function  (PDF)  given  by  Eq.  (5)  in  Section  2.3. 

A  model  developed  from  this  study  for  Sy  caused  by  scattering  from  suspended 
glacial  silt  particles  (and  thus  applicable  to  the  Kenai  River)  is  presented  in  the 
appendix.  The  predictions  of  the  model  agree  well  with  our  estimates  of  Sy,  but 
the  effect  of  possible  additional  scattering  firom  bubbles  needs  to  be  determined. 
This  model  should,  however,  motivate  further  study  on  the  relationship  between  Sy 
and  turbidity,  and  it  may  be  particularly  useful  for  acoustical-limnological  studies  in 
Alaskan  lakes  containing  glacial  meltwater. 
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Our  study  quantifies  how  using  horizontal-looking  sonars  to  detect  fish  in  rivers 
is  reverberation  limited  and  not  noise  limited.  Knowing  Sy,  for  example,  we  can  then 
estimate  the  minimum  target  strength  of  fish,  T^mm,  that  can  be  reliably  detected. 
Assuming,  for  example,  that  T<Smin  niust  be  at  least  10  dB  more  than  the  equivalent 
target  strength  of  the  volume  of  water  ensonfied,  TSyou  fhen  it  follows  that  TSmm  > 
TSyoi  -f  10  dB,  where 

rSv.l  =  101ogj„(yi?"V’)+S.. 

This  equation  can  be  used  to  evaluate  the  impact  of  changing  the  pulse  length,  r,  of 
the  sonar  beam  as  parameterized  by  ip,  and  range,  R,  with  respect  to  fish-detection 
goals  and  algorithms. 

Another  important  estimate  that  depends  on  Sy  is  the  root  mean  square  (rms) 
error  in  angular  location,  ^rjns,  when  tracking  fish  with  split-beam  sonars.  This  er¬ 
ror  is  proportional  to  where  SNR  is  the  signal-to-reverberation  ratio  and  is 

determined  by 

where  TS^h  is  the  nominal  target  strength  of  a  fish  at  side  aspect. 
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1  INTRODUCTION 


In  several  important  commercial,  and  sport,  fishing  salmon  river  systems,  the 
Alaska  Department  of  Fish  and  Game  (ADF&G)  uses  side-looking  sonars  mounted 
along  riverbanks  to  count  adult  salmon  during  their  upstream  migration.  During  the 
1990s,  the  average  commercial  catch  of  Pacific  salmon  in  Alaska  was  175  million  fish; 
the  commercial  catch  exceeded  215  million  fish  twice  in  that  decade.  High  harvest 
levels  of  Pacific  salmon  from  the  mid  1970s  to  the  present  have  been  maintained  by 
a  combination  of  careful,  active  management  [1]  and  favorable  environmental  condi¬ 
tions  [2], [3].  In  most  cases,  the  harvests  have  been  managed  to  achieve  approximately 
fixed  spawning  stock  sizes  that  do  not  change  with  year-to-year  variation  in  abun¬ 
dance.  A  major  technical  challenge  has  been  to  estimate  spawning  stock  size  while 
the  runs  and  the  fisheries  are  actively  under  way.  Some  runs  last  just  a  few  weeks  and 
involve  tens  of  millions  of  adult  salmon  moving  into  vast,  remote  wilderness  areas  to 
spawn.  The  fisheries  managers  can  visually  estimate  or  monitor  passage  with  aerial 
observations,  counting  towers,  or  counts  through  weirs  for  systems  with  clear  water. 
Sonar  has  been  used  since  the  1960s  for  important  spawning  systems  in  Alaska  that 
are  glacially  influenced  or  otherwise  not  suitable  for  visual  observation  [4], [5].  Early 
studies  involved  arrays  of  up-looking  sonar  produced  by  the  Bendix  Corporation.  By 
1965  engineers  from  Bendix  tested  a  side-looking  design  that  formed  the  basis  of  the 
machines  that  are  still  used  in  many  important  salmon-producing  rivers  in  Alaska.  In 
1997,  the  Alaska  Department  of  Fish  and  Game  decided  to  replace  the  1960s  vintage 
sonars  with  split-beam  systems.  Field  tests  of  side-looking  split-beam  sonar  began  in 
1998  and  continued  in  1999. 

This  report  summarizes  results  of  a  project  involving  the  Applied  Physics  Lab¬ 
oratory  (APL),  ADF&G,  Biosonics  Inc.,  and  Hydroacoustics  Technology  Inc.  (HTI) 
to  measure,  and  model,  the  background  acoustic  reverberation  in  two  Alaskan  rivers. 
Background  reverberation  is  defined  as  the  average  acoustic  intensity  of  the  echo 
versus  time  (range)  in  the  absence  of  any  contribution  from  fish.  The  rivers  were 
the  Kenai  River,  located  on  the  Kenai  Peninsula,  and  the  Wood  River,  which  drains 
into  Bristol  Bay  near  Dilhngham.  The  measurements  were  made  at  two  frequencies: 
420  kHz,  using  a  system  built  by  Biosonics  Inc.,  and  200  kHz,  using  a  system  built  by 
HTI.  For  each  system,  a  combination  of  transducers  and  pulse  lengths  was  employed 
as  discussed  in  the  following  sections. 

Background  reverberation  determines,  in  part,  the  minimum  size  class  of  fish  that 
is  detectable  with  sonar  and  the  accuracy  of  fish-tracking  and  counting  algorithms  that 
use  the  complex  phase  of  the  fish  echo,  and  is  necessary  to  establish  rigorous  detection 
probabilities  for  all  fish  size  classes.  Furthermore,  with  a  reasonably  accurate  model 
for  background  reverberation,  it  is  possible  to  simulate  the  acoustic  background  to 
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evaluate  the  consequences  of  a  change  in  sonar  system  parameters  such  as  pulse  length, 
transducer  beamwidth,  or  frequency  on  fish-detection  and  tracking  algorithms.  We 
determined  that  the  background  reverberation  in  both  riverine  environments  studied 
was  set  by,  in  order  of  importance,  (1)  scattering  from  bubbles  and  particles  of  glacial 
silt  entrained  in  the  water  column  {volume  backs cattering),  and  (2)  scattering  from 
the  rough  surface  of  the  riverbed  {bottom  backscattering). 

The  report  is  organized  as  follows.  Section  2  summarizes  the  measurements  made 
at  the  Kenai  site  and  the  resulting  sonar  models  for  background  reverberation,  and 
Section  3  repeats  this  summary  for  the  Wood  River  site.  Note  that  Section  2  is 
considerably  longer  than  Section  3  as  it  includes  introductory  remarks  and  a  discussion 
of  important  notation  used  throughout  the  report.  In  Section  4  we  summarize  and 
present  conclusions,  and  in  the  appendix  we  present  a  model  for  volume  scattering 
from  glacial  silt  that  applies  to  the  Kenai  River. 
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2  MEASUREMENTS  IN  THE  KENAI  RIVER 

2.1  Site  Description  and  Geometry 

The  Kenai  River  measurements  were  made  8  and  9  June  1999  at  a  site  8.5  miles 
upstream  from  the  mouth  of  the  river  (Fig.  1)  on  the  right  bank  that  is  used  by 
ADF&G  for  sonar  counting  of  migrating  chinook  salmon  {Oncorhyncus  tshawytscha). 
The  Kenai  River  at  this  site  is  approximately  90  m  wide  and  5  m  deep  at  the  thalweg 
at  mean  tide  levelj  the  river  depth  undergoes  a  tidal  fluctuation  of  approximately 
±1.5  m.  Additional  details  concerning  this  site  are  available  in  Refs.  6  and  7.  The 
HTI  system  was  located  about  50  m  upstream  from  the  ADF&G  sonar  counting  site, 
and  the  Biosonics  system  was  located  another  10  m  upstream  from  the  HTI  sys¬ 
tem.  Fathometer  readings  were  made  from  cross-river  transects  at  both  locations  the 
day  before  measurements  started.  These  data  were  used  to  construct  a  bathymetric 
profile  of  the  right  bank  of  the  Kenai  River,  starting  from  each  transducer  location. 
Figure  2,  for  example,  shows  the  position  of  the  Biosonics  transducer  along  with 
important  geometric  variables,  which  are  the  transducer’s  height  above  the  bottom, 
Ht  =  0.48  m,  the  average  bottom  slope,  9s,  and  the  water  depth,  D,  versus  range, 
R,  from  the  sonar.  At  the  Biosonics  location  the  bottom  slope  varied  from  about 
4°  just  out  from  the  transducer  to  about  3°  at  a  range  of  20  m,  with  3.5°  represent¬ 
ing  a  nominal  average  slope  of  the  riverbed.  The  fathometer  readings  made  at  the  HTI 


Figure  1.  Map  of  the  lower  Kenai  River  showing  location  of  the  sonar  test  site. 
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Range  (m) 


Figure  2.  Experimental  geometry  for  acoustic  reverberation  measurements  made 
in  the  Kenai  River  and  key  geometric  variables.  View  represents  a  portion  of  the 
river  cross  section  at  its  right  bank  extending  out  to  the  thalweg  corresponding  to 
the  location  of  the  Biosonics  system.  The  transducer’s  position  is  shown  by  the  large 
triangle  at  range  =  0  m.  The  bottom  contour  has  been  constructed  from  fathometer 
readings  made  during  cross-river  transects. 

location  showed  a  similar  bathymetry,  but  the  average  slope  across  the  transect  was 
closer  to  3°;  also,  for  the  HTI  system  Hr  —  0.25  m.  Note  that  D  changed  over  the 
time  period  (approximately  3  hours)  during  which  the  acoustic  measurements  were 
made  owing  to  tidal  effects.  We  obtained  the  correct  D  at  the  time  of  each  acoustic 
measurement  from  depth  measurements  made  by  ADF&G  personnel  at  their  sonar 
site. 

Because  the  Kenai  River  undergoes  tidal  fluctuations  at  the  experimental  site,  it 
was  important  to  first  evaluate  the  effect,  if  any,  of  these  fluctuations  on  the  acoustic 
environment.  Specifically,  it  is  possible  that  an  incoming  tidal  flow  could  transport 
with  it  high-salinity  seawater  that  would  mix  with  the  less  saline  river  water,  increas¬ 
ing  both  the  acoustic  attenuation  and  the  sound  speed.  To  evaluate  this,  several 
conductivity-temperature-depth  (CTD)  measurements  were  made  with  a  Sea-Bird 
CTD  unit  on  7  and  8  June  to  sample  various  tidal  phases.  Figure  3  shows  the  results 
of  a  CTD  measurement  made  at  1030  on  8  June,  which  is  representative  of  the  other 
measurements.  (Note:  All  times  reported  here  are  local.)  The  three  types  of  lines 
represent  CTD  deployments  made  at  three  points  roughly  evenly  spaced  across  the 
river.  The  dotted-dashed  line  corresponds  to  readings  closer  to  the  right  bank;  the 
solid  line  corresponds  to  readings  near  the  thalweg,  and  therefore  extends  the  deepest; 
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the  dashed  line  corresponds  to  readings  closer  to  the  left  bank.  The  measurements 
show  that  (1)  salinity  was  extremely  low  (approximating  fresh  water),  and  (2)  the 
sound  speed  was  primarily  determined  by  water  temperature,  which  was  nearly  con¬ 
stant  with  depth.  During  this  particular  measurement  the  water  temperatme  was 
about  8.25°C,  but  for  others  the  temperature  ranged  between  8°  and  9°C,  which 
puts  the  sound  speed,  c,  in  the  range  1438.5  m/s  to  1442.6  m/s  [8].  No  salinity  of 
any  acoustic  significance  was  measured,  and  thus  the  main  effect  of  the  tide  on  the 
acoustic  environment  was  to  vary  the  water  depth,  D. 


Sound  Speed  (m/s)  Temperature  (deg  C)  Salinity  (PSU) 

Figure  3.  Results  of  the  CTD  measurements  made  in  the  Kenai  River  at  1030  local 
time  on  8  June  1999.  Dotted-dashed  line  is  the  measurement  made  approximately 
20  m  out  from  the  right  bank,  the  solid  line  is  the  measurement  made  near  the 
thalweg,  and  the  dashed  line  is  the  measurement  made  approximately  20  m  out 
from  the  left  bank.  The  salinity  is  displayed  here  in  practical  salinity  units,  or  PSU. 


2.2  Acoustic  Measurements 

One  of  the  key  objectives  of  the  1999  field  effort  on  the  Kenai  and  Wood  rivers  was 
to  study  sonar  fish-tracking  and  enumeration  strategies  using  the  Biosonics  and  HTI 
systems.  This  report  focuses  only  on  the  environmental  acoustic  measurements  made 
with  these  two  systems  in  order  to  understand  background  reverberation  in  sonar- 
based  fish  counting.  However,  a  protocol  had  to  be  designed  to  measure  acoustic 
reverberation  in  a  convenient  manner  while  also  utilizing  the  optimum  positioning  of 
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the  two  sonars  for  fish  detection.  This  protocol,  which  was  also  followed  during  the 
Wood  River  measurements,  is  outlined  below. 

Each  operator  determined  the  optimum  transducer  pitch  angle,  6p  (defined  as  pos¬ 
itive  downward),  for  counting  fish  at  their  particular  site.  Typically,  the  transducers 
need  to  be  aimed  downward  a  few  degrees  from  horizontal  such  that  the  main  lobe  of 
the  transducer  beam  is  placed  within  the  path  of  the  migrating  salmon.  The  salmon 
swim  within  a  few  centimeters  of  the  riverbed  to  take  advantage  of  the  reduced  water 
velocity  within  the  bottom  boundary  layer  [9].  The  final  established  pitch  angle,  call 
it  Op  =  Opp  for  fish  angle,  is  influenced  by  the  transducer’s  beamwidth,  the  bottom 
slope,  and  the  composition  and  granularity  of  the  riverbed.  The  operator  selects  Opp 
by  increasing  Op  imtil  the  ensuing  backscatter  from  the  riverbed  increases  significantly 
and  begins  to  dominate  the  background  reverberation.  At  this  point  Op  is  reduced 
perhaps  a  degree  or  two  (i.e.,  back  toward  the  horizontal)  with  subsequent  reduction 
in  the  background  reverberation.  Once  Opp  was  established,  background  acoustic  re¬ 
verberation  measurements  were  made  at  this  angle  and  at  various  pitch  angles  above 
and  below  Opp.  For  the  Biosonics  system  Opp  was  determined  to  be  4°;  the  Op  settings 
used  for  the  acoustic  measurements  included  this  angle  as  well  as  —2°,  0°,  5°,  and 
8°.  For  the  HTI  system  Opp  was  determined  to  be  2°,  and  the  Op  settings  included 
this  angle,  0°,  and  4°.  (The  transducer  pitch  angles  reported  here  are  to  the  nearest 
degree.) 

A  measurement  was  defined  as  100  pings  transmitted  at  1-s  intervals.  Upon  com¬ 
pletion  of  a  measurement,  the  data  were  examined  briefly,  and  then  the  measurement 
was  repeated  using  a  different  Op  or  pulse  length,  r.  For  the  420-kHz  system,  measme- 
ments  were  made  at  r  =  0.2  ms  and  0.4  ms.  For  the  200-kHz  system,  measurements 
were  made  at  r  =  0.2  ms  and  0.5  ms  and  with  an  FM  pulse  with  a  bandwidth  of 
10  kHz  and  an  effective  pulse  length  of  0.18  ms.  Other  key  sonar  system  parameters 
used  during  the  Kenai  measurements  (which  are  discussed  in  the  text)  are  summarized 
in  Table  1. 

Table  1.  Sonar  system  parameters  used  for  acoustic  measurements 

in  the  Kenai  River. 


System 

Frequency 

Source  Level 

(hw 

4>n 

10  logic  V’ 

* 

(kHz) 

(dB) 

(deg) 

(deg) 

(dB) 

Biosonics  (circ.) 

420 

216.1 

2.3 

2.3 

-30.4 

HTI  (ellip.) 

200 

217.6 

10.1 

2.3 

-24.0 
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Reverberation  Due  to  Volume  Backscattering 


Volume  reverberation  is  the  term  used  to  describe  the  scattering  of  sound  back 
toward  the  transducer  from  the  total  volume  of  water  ensonified  by  the  sonar  pulse. 
Here  we  introduce  important  notation  used  in  both  the  analysis  and  the  reporting  of 
results.  The  discussion  follows  closely  the  one  given  by  Dahl  [8]. 

Consider  a  cloud  of  scatterers  at  range  R  corresponding  to  the  cloud’s  center.  If 
we  ignore  for  now  the  loss  due  to  sound  attenuation  in  water,  an  elemental  volume, 
dV,  produces  a  backscattered  intensity  at  the  receiver  of  dibs,  given  by 


dibs  — 


linc^v^^ 


The  quantity  SydV  assumes  the  role  of  the  backscattering  cross  section,  abs,  for  an 
assemblage  of  scatters  within  elemental  volume  dV,  where  Sy  is  the  backscattering 
cross  section  per  cubic  meter  of  water  and  has  units  of  m“^.  The  total  backscattered 
intensity  results  from  summing  all  dV,  some  of  which  are  away  from  the  acoustic 
axis,  to  form  the  effective  scattering  volume,  V.  For  these  off-axis  contributions,  the 
incident  and  backscattered  intensity  axe  reduced  slightly  according  to  the  intensity 
beam  pattern  b{d,4>).  The  net  effect  leads  to  the  concept  of  an  effective  volume, 
or  reverberation  volume  [10],  based  on  integration  of  the  two-way  intensity  pattern 
6^(0,  (f).  If  we  define  ip  as  the  integral  of  5^(0,  (j))  over  all  solid  angles,  then  the  effective 


scattering  volume  at  range  R  for  a  pulse  of  length  r  is  V 
backscattered  intensity  is 

r  ^orl  „  cr  _2 


^Rrtl),  and  the  total 


Ihs  — 


where  the  incident  intensity  is  referenced  back  to  Iq  via  spherical  spreading  with 
line  =  loi'f'o/R)^  and  Tq  is  a  reference  distance  set  to  1  m.  A  convenient  expression 
[11]  for  ^  which  we  use  in  this  analysis  is 

where  and  (pN  are  the  full,  one-way  —3  dB  beamwidths  (expressed  in  radians) 
for  the  wide  (W)  and  narrow  (N)  angles  that  describe  the  elliptical  HTI  transducer 
beams  used  in  this  study.  We  also  use  Eq.  (3)  for  the  circular  Biosonics  transducer 
beams  by  setting  (f>w  =  <I>n-  Values  oi^w,  <I>n,  (in  degrees)  and  rp  for  the  two  systems 
during  the  Kenai  measurements  are  given  in  Table  1. 

The  sonar  equation  for  volume  reverberation  is  the  decibel  equivalent  of  Eq.  (2): 

RLy  =  SL- 40  logic  R  ~  -I-  S'v  -f  10  logic 
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where  SL  is  the  source  level,  RL^  is  the  reverberation  level  specifically  attributed  to 
volume  scattering  and  is  equal  to  lOlog^o/fts)  ^-nd  the  effect  of  two-way  absorption 
loss  is  included  as  2aR,  where  a  is  the  water  attenuation  coefficient  (in  dB/m).  Given 
an  effective  scattering  volume  V,  lOlogio(svl^)  is  the  equivalent  target  strength  (TS) 
of  the  scattering  volume,  which  is  comparable  directly  to  the  TS  of  a  fish.  The 
volume  scattering  strength,  Sv,  is  10  log  (in  dB  re  1  m“^);  this  is  the  key  variable 
pertaining  to  voliime  scattering  that  we  shall  estimate  and  report.  Note  that  Sy  can 
vary  with  both  range  and  depth.  However,  with  the  exception  of  episodic  boat  wakes, 
a  range-  and  depth-independent  Sy  is  a  suitable  approximation  for  the  two  riverine 
environments  studied  here. 

Note  that  the  two  systems  were  necessarily  configured  differently  in  order  to  make 
the  kind  of  measurements  demanded  by  this  study.  However,  with  each  we  ultimately 
obtained  a  100-ping  ensemble  average  of  calibrated  acoustic  intensity  for  a  given  time- 
delay  bin;  this  was  equated  to  range  using  a  nominal  sound  speed  of  1440  m/s.  In 
particular,  we  worked  with  ensemble-averaged  intensity  from  the  200-kHz  system, 
whereas  we  worked  with  the  entire  set  of  100  pings  from  the  420-kHz  system  and  com¬ 
puted  ensemble  averages  during  postprocessing.  The  different  configurations  produce 
a  qualitative  difference  in  the  data  display,  e.g.,  some  small  range  averaging  is  already 
included  in  the  200-kHz  data,  but  otherwise  these  differences  are  of  no  consequence. 
The  exception  is  when  an  occasional  fish  was  detected  in  some  of  the  measurements, 
and  could  not  be  removed  from  the  200-kHz  ensemble-averaged  data.  In  any  case,  the 
primary  observable  from  these  measurements  is  the  ensemble-averaged  backscattered 
acoustic  intensity  versus  echo  time  delay.  We  report  this  quantity  as  the  mean  rever¬ 
beration  level,  or  RL  without  any  subscript,  when  expressed  in  dB  with  reference  to 
1  //Pa.  (The  shorthand  notation  dB  re  1  //Pa  is  hereafter  used.)  Fluctuations  within 
the  100-ping  ensemble  were  also  studied  using  data  from  the  420-kHz  system. 

To  estimate  Sy  we  first  computed  RL  using  data  from  measurements  where  the 
transducer  orientation  was  close  to  horizontal.  At  this  orientation,  the  sonar  beam 
ensonifies  only  a  portion  of  the  water  column  in  a  given  range  bin,  and  backscatter 
from  the  surface  and  bottom  boundaries  is  negligible.  We  thus  assume  that  the 
measured  RL  is  due  solely  to  volume  reverberation,  and  obtain  an  estimate  of  Sy 
by  solving  for  this  quantity  in  Eq.  (4).  We  remark  that  Sy  then  is  also  the  result  of 
linearly  averaging  acoustic  intensity  and  then  taking  the  decibel  equivalent. 

Our  estimation  procedure  is  illustrated  in  the  following  two  figures  showing  mea¬ 
surements  made  at  420  kHz  with  Op  set  to  0°  (Figure  4)  and  at  200  kHz  with  Op  set 
to  2°  (Figure  5).  Both  measurements  were  made  at  approximately  noon  on  8  June. 
Figure  4  (top)  is  an  image  of  the  100  pings  versus  ping  number,  or  time  in  seconds, 
and  range  from  the  sonar.  The  image  shows  a  dramatic  increase  in  RL  at  about 
ping  23  which  persists  for  more  than  80  s.  This  is  caused  by  the  bubbly  wake  of  a 
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Figure  4.  Top:  Image  of  100  pings  transmitted  at  1-s  intervals  and  a  frequency 
of  420  kHz  versus  ping  number  and  range  from  sonar,  with  lighter  shade  of  gray 
indicating  higher  reverberation  level.  The  wake  of  a  passing  motor  boat  can  be  seen 
starting  at  ping  23,  and  the  cross-river  dimension  of  the  wake  gradually  increases 
from  about  2  m  to  10  m  in  40  s.  Bottom:  Reverberation  level  {RL)  measured  during 
ping  number  33  (thin,  erratic  line),  mean  RL  based  on  the  100-ping  ensemble  average 
(thick,  gray  line),  and  a  model  for  RL  (thick,  black  line). 

passing  boat  which  was  also  noted  in  our  field  observations.  The  wake  broadening 
caused  by  turbulent  diffusion  can  be  observed,  and  the  cross-river  dimension  of  the 
wake  increases  from  about  2  m  to  10  m  in  40  s.  Figure  4  (bottom)  shows  the  calibrated 
reverberation  level  during  ping  number  33  (thin,  erratic  line),  plotted  along  with  the 
mean  RL  (thick,  gray  line)  based  on  the  100-ping  ensemble  average,  and  a  model 
for  RL  (thick,  black  line).  The  model  is  Eq.  (4)  evaluated  using  the  source  level, 
a  =  0.060  dB/m  [8],  and  i/)  values  from  Table  1  and  with  Sy  set  to  -58  dB,  our 
estimate  of  Sy  at  420  kHz  for  this  particular  run.  If  we  ignore  the  episodic  wake 
feature,  which  increases  Sy  by  about  15  dB,  our  estimate  of  Sy  is  clearly  supported 
by  the  data. 

Figure  5  shows  the  equivalent  ensemble-averaged  RL  versus  range  measured  with 
the  200-kHz  sonar.  Here,  two  runs  taken  within  4  min  of  each  other  are  shown  by  the 
two  gray  lines  that  appear  as  one,  and  a  third  run  taken  about  10  min  later  is  shown 
by  the  thin  black  line.  A  model  for  the  mean  RL  (thick,  black  line)  again  represents 
an  evaluation  of  the  sonar  equation  using  the  source  level,  a  =  0.013  dB/m  [8],  and 
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Figure  5.  Mean  RL  at  200  kHz  versus  range  from  sonar  based  on  a  100-ping 
ensemble  average  (thick  gray  line)  and  a  model  for  RL  (thick  black  line).  The  thin 
black  line  is  another  measurement  of  RL  taken  approximately  10  min  later. 


Ip  values  from  Table  1;  however,  here  we  have  set  Sy  to  —71  dB,  our  estimate  of  Sy  at 
200  kHz.  The  jump  in  measured  RL  near  a  range  of  20  m  is  likely  due  to  an  unknown 
scattering  feature  such  as  a  piece  of  stationary  debris;  otherwise  the  estimate  of  Sy 
is  well  supported  by  the  data.  An  apparent  exception  to  this  agreement  begins  at  a 
range  of  25  m.  This  range  corresponds  to  the  bubble-front  of  the  boat  wake  seen  in 
Figure  4.  There  is  some  evidence  from  studies  of  ship  wakes  that  sound  absorption 
within  wakes  is  more  severe  at  200  kHz  than  at  420  kHz.  However,  we  have  no 
knowledge  of  the  exact  position  of  the  two  sonar  beams  with  respect  to  the  wake 
structure.  The  420-kHz  beam,  for  example,  could  have  been  ensonifying  the  more 
diflFuse  edge  of  the  wake,  while  the  200-kHz  beam  could  have  been  directed  toward 
the  more  bubble-dense  interior  of  the  wake.  In  any  case,  a  repeat  of  the  200-kHz 
measurement  approximately  10  min  later  showed  the  RL  in  the  range  20-30  m  had 
returned  closer  to  the  values  predicted  by  our  model  using  Sy  =  —71  dB. 

Notwithstanding  the  differences  seen  in  the  response  to  episodic  boat  wakes,  om 
estimates  of  the  Sy  at  420  kHz  are  about  13  dB  greater  than  those  at  200  kHz,  which 
implies  that  Rayleigh  scattering  is  in  effect,  where  the  scattered  intensity  goes  as 
the  fourth  power  of  frequency.  In  the  appendix  we  present  a  model  for  Rayleigh 
scattering  from  suspended  particles  associated  with  the  glacial  silt  load  carried  in  the 
Kenai  River. 

Figure  6  shows  a  sample  histogram  of  RL  at  a  range  of  15  m,  computed  for 
each  ping  of  the  100-ping  420-kHz  measurement  shown  in  Figure  4.  The  curve  is  a 
probabihty  density  function  (PDF)  and  is  given  by  [12] 

=  —i^x — '  (5) 
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Figure  6.  Sample  histogram  of  the  ping-by-ping  420-kHz  reverberation  level  (in 
dB)  corresponding  to  a  range  of  15  m  in  the  data  shown  in  Figure  4.  Curve  is  a 
model  PDF,  and  significance  probability  of  the  statistic  is  greater  than  0.75. 

where  the  random  variable  s  represents  acoustic  reverberation  level,  ki  =  lOlogioC, 
and  A  is  the  mean  acoustic  reverberation.  A  chi-square  (x^)  goodness-of-fit  test 
comparing  the  model  PDF  with  the  data  shown  in  Figure  6  results  in  a  significance 
probability  of  the  x^  statistic  being  greater  than  0.75,  suggesting  the  PDF  is  a  very 
good  model  for  the  observed  fluctuations.  Since  for  a  given  range,  RL  and  are 
related  through  a  constant  offset,  we  can  thus  say  that  fluctuations  in  Sv  for  a  given 
range  are  also  distributed  according  to  the  PDF  given  by  Eq.  (5). 

The  measurements  shown  in  Figures  4  and  5  each  took  100  s  to  make,  and  as¬ 
sumptions  regarding  the  stationarity  of  in  order  to  compare  measurements  with 
the  model  PDF  of  Eq.  (5)  are  easily  satisfied.  Over  the  3-hour  duration  of  the  ex¬ 
periment,  however,  we  observed  a  significant  variation  in  our  estimates  of  Sy  The  30 
measurements  of  Sv  made  with  the  200-kHz  system,  for  example,  were  roughly  evenly 
distributed  between  a  high  of  *-62  dB  and  low  of  -75  dB.  To  obtain  a  measure  of  the 
central  tendency  of  S^,  we  choose  to  average  these  estimates  using  101ogio((10^’'/^°)), 
with  the  result  being  —  67  dB.  This  same  average  performed  on  the  estimates  made  at 
420  kHz  yields  —53  dB.  The  difference,  about  14  dB,  remains  close  to  that  predicted 
by  Rayleigh  scattering.  Finally,  we  use  ±  4  dB  as  a  measure  of  the  expected  spread 
in  Sy  for  both  the  200-kHz  and  420-kHz  measurements;  this  range  encompasses  22 
of  the  30  measurements  made  at  200  kHz,  and  11  of  the  15  measurements  made  at 
420  kHz. 
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2.4  Reverberation  Due  to  Bottom  Backscattering 

Bottom  reverberation  is  the  term  used  to  describe  aooustic  scattering  in  the 
direction  back  toward  the  transducer  from  the  surface  of  the  riverbed  ensonified  by 
the  sonar  pulse.  Some  additional  notation  is  now  introduced,  starting  with  the  sonar 
equation  for  backscattering  from  the  riverbed: 

RLb  =  SL- 40  log  R  -  2aR  +  56  +  10  logjo  A- B.  (6) 

In  Eq.  (6),  A  is  the  effective  area  of  the  riverbed  contributing  to  the  backscatter,  RL^ 
is  the  reverberation  level  specifically  associated  with  riverbed  scattering,  and  B  is  the 
round-trip  beam  loss  (in  dB).  Given  an  effective  area.  A,  then  lOlogio  A  +  Sb  is  the 
equivalent  TS  of  an  area  of  riverbed  with  scattering  strength  Sb  (in  dB),  and  it  is  the 
key  variable  pertaining  to  bottom  scattering  that  we  shall  estimate  and  report. 

Note  that,  unfike  for  5„,  estimating  Sb  using  sonar  beams  that  are  oriented  hori¬ 
zontally  is  a  delicate  task.  The  reverberation  caused  by  scattering  from  the  riverbed 
reaches  the  transducer  through  angles  away  from  that  of  the  transducer’s  main  lobe, 
and  the  scattering  from  the  riverbed  usually  originates  at  very  small  grazing  angles, 
9a-  These  factors  tend  to  produce  a  magnified  response  to  very  small  changes  in  the 
scattering  geometry,  making  it  more  difficult  to  obtain  robust  estimates  of  Sb-  This 
situation  is  improved  somewhat  by  pitching  the  transducer  downward,  at  angle  9p, 
so  that  more  of  the  main  lobe  is  pointed  toward  the  riverbed.  Figure  7  shows  this 
configuration  using  a  simplified  representation  of  a  transducer  main  lobe.  A  patch 
on  the  riverbed  (marked  A)  scatters  sound  back  to  the  transducer  at  grazing  an¬ 
gle  Oq.  For  a  given  composition  and  roughness,  the  scattering  strength,  Sb,  of  the 
riverbed  is  a  function  only  of  9c,  which  for  the  geometry  shown  in  Figure  7  is  set  by 

sin  9a  —  Hp  sin(90  +  9s)  /  R-  (7) 

The  RLb  recorded  by  the  transducer,  however,  depends  on  the  effective  launch  angle, 
9'j^,  linked  to  the  scattering  patch.  For  geometric  considerations,  it  is  convenient  to 
define  a  vertical  launch  angle,  9^,  with  respect  to  the  horizontal,  and  to  define  9'^^  with 
respect  to  the  transducer’s  acoustic  axis,  which  puts  9l  —  9a  +  9$  and  9'^  =  9l  —  9p. 
The  round-trip  beam  loss  is  given  hy  B  =  — lOlogio^^(^L)O))  where  h'^{9,(j))  is  the 
two-way  intensity  beam  pattern  as  a  function  of  vertical  {9)  and  azimuthal  (0)  angles, 
,and  the  effective  area  A  is  given  by 

To  evaluate  B,  and  the  integral  in  Eq.  (8),  we  construct  numerical  representations  of 
using 
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Figure  7.  The  geometry  of  backscattering  from  a  sloping  riverbed.  Scattering  orig¬ 
inates  from  a  patch  of  riverbed  (thick  black  line  labeled  A)  at  the  grazing  angle  0g- 
The  slope  of  the  riverbed,  9s,  the  transducer  pitch  angle,  9p,  and  the  transducer’s 
height  above  the  bottom.  Hr,  are  the  key  variables  that  determine  the  launch  angle, 

Ol,  and  the  effective  launch  angle,  9'j^,  and  thereby  the  reverberation  level  recorded 
by  the  transducer. 

b‘^{9,4>)  ^  [bN{0)bw{<f>)]‘^-  (9) 

This  approximation  is  based  on  the  assumption  that  a  one-way  intensity  pattern  de¬ 
fined  in  both  the  9  and  0  planes  can  be  specified  by  multiplying  two  single-coordinate 
one-way  intensity  patterns  [13].  The  necessary  one-way  intensity  patterns  for  bN{9) 
and  bw{<f>)  were  measured  by  Biosonics  and  HTI  for  the  transducers  used  in  this  study, 
and  the  angular  resolution  in  these  measurements  was  increased  by  interpolation  for 
the  purpose  of  carrying  out  the  integral  in  Eq.  (8). 

Figure  8  shows  the  mean  RL  for  six  sets  of  measurements  at  200  kHz;  again,  each 
set  represents  the  ensemble  average  of  100  pings  for  which  9p  =  4°.  Each  measurement 
set  was  separated  by  2  min;  two  sets  were  made  using  a  cw  transmission  length  of 
r  =  0.2  ms,  two  sets  were  made  using  a  cw  transmission  length  of  r  =  0.5  ms,  and 
two  sets  made  using  an  FM  pulse  with  an  effective  r  of  0.18  ms.  The  data  have  been 
normalized  to  an  equivalent  r  of  0.2  ms  by  subtracting  3.98  dB  from  the  measurements 
made  with  r  =  0.5  ms  and  adding  0.46  dB  to  the  measurements  made  with  an  effective 
r  of  0.18  ms.  The  six  measurements  agree  well  following  this  normalization,  which 
shows  that  RL  increases  or  decreases  according  to  lOlog^o  y  (in  dB).  The  dotted- 
dashed  line  is  modeled  RL^  when  using  =  -71  dB,  and  the  dashed  line  is  modeled 
RLf,  when  using  a  scattering  strength  model  (in  dB)  given  by 
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Figure  8.  Mean  RL  measured  at  200  kHz  versus  range  from  sonar  based  on  100- 
ping  ensemble  average  (gray  line);  model  for  RLy  (dotted-dashed  line);  model  for 
RLh  (dashed  line),  and  the  incoherent  sum  of  RLy  and  RLi,  (solid  line). 

'S'6(0g)  =  401ogio(sm6>G) +  2.  (10) 

The  solid  black  line  is  the  incoherent  sum  of  these  two  components  as  computed  by 

101ogio(10^^’"^^°+^^'’/^°).  (11) 

Both  the  models  and  the  data  show  that  when  the  transducer  is  pitched  downward  to 
dp  =  4°  scattering  from  the  riverbed  dominates  RL  between  ranges  of  roughly  5  and 
10  m,  whereas  volume  scattering  tends  to  dominate  outside  this  region.  The  models 
result  from  evaluating  the  two  sonar  equations,  Eq.  (4)  and  Eq.  (6),  using  data  from 
Table  1  and  the  applicable  beam-pattern  functions  6jv(^)  and  bw{(l>),  with  Op  —  4° 
and  Hp  =  0.25  m  being  key  inputs  for  the  evaluation  of  RLb. 

Figure  9  shows  the  mean  RL  for  two  sets  of  measurements  made  at  420  kHz  for 
which  Op  =  8°.  Each  measurement  set  was  also  separated  by  2  min.  One  set  was 
made  using  a  cw  transmission  length  of  r  =  0.2  ms,  and  the  other  set  with  r  =  0.4  ms, 
which  has  been  normalized  to  an  equivalent  r  =  0.2  ms  by  subtracting  3  dB.  The 
dotted-dashed  line  is  the  modeled  RLy  when  using  Sy  =  —53  dB,  and  the  dashed  line 
is  the  modeled  RLb  when  using  a  scattering  strength  model  (in  dB)  given  by 

SbiOa)  =  40 logio(sin Oq) +7  .  (12) 

The  sohd  black  line  represents  the  incoherent  sum  of  these  two  components.  As  in 
the  example  shown  in  Figure  8,  Op  =  8°  and  Hp  =  0.48  m  are  important  inputs  for 
the  evaluation  of  RLb-  In  this  case  both  the  models  and  the  data  show  that  when  the 
transducer  is  pitched  downward  to  Op  =  8°  scattering  from  the  riverbed  dominates  the 
RL  between  ranges  of  roughly  5  and  20  m,  whereas  volume  scattering  tends  to  dom¬ 
inate  outside  this  region.  Note  that  the  data  between  5  and  20  m,  and  also  the  data 
shown  in  Figure  8  between  5  and  10  m,  display  a  slight  oscillatory  character  because 
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Figure  9.  Mean  RL  measured  at  420  kHz  versus  range  from  sonar  based  on  100- 
ping  ensemble  average  (gray  line);  model  for  RL^  (dotted-dashed  line);  model  for 
RLb  (dashed  line),  and  the  incoherent  sum  of  RLv  and  RLb  (solid  line). 

the  mean  RL  that  originates  from  the  riverbed  does  not  represent  a  true  ensemble, 
since  the  same  patch  of  riverbed  is  ensonified  with  each  ping.  Nevertheless,  given  this 
deficiency  in  the  measurement,  our  models  for  RL  still  tend  to  agree  approximately 
with  the  data  in  each  case. 

The  two  empirical  models  for  Sb  were  derived  first  by  examining  estimates  of  Sb 
obtained  by  inverting  Eq.  (6),  and  then  mapping  the  slant  range,  R,  to  the  grazing 
angle  6g  to  get  Sb  versus  9g-  We  chose  to  model  Sb  using  a  Bragg  scattering  law, 
where  SbiOa)  ~  4OlogioSin(0G)-  However,  only  a  limited  range  of  grazing  angles  is 
involved  in  scattering  from  the  riverbed.  For  example,  the  region  of  definitive  riverbed 
scattering  shown  in  Figure  9  between  about  5  m  and  20  m  involves  grazing  angles 
between  only  2°  and  4°.  It  is  quite  possible  that  other  candidate  scattering  laws 
could  apply,  e.g.,  Lambert’s  law,  where  SbiOc)  ~  20 logic  sin(0G).  When  comparing 
the  resultant  total  RL  with  the  measured  RL,  though,  the  40  logic  sin(^G)  dependence 
produced  the  most  acceptable  result.  Were  Bragg  scattering  to  be  the  true  modality 
for  acoustic  backscattering  from  the  riverbed,  then  this  scattering  would  originate 
from  roughness  at  very  small  scales  corresponding  to  the  Bragg  wavelength,  which 
is  half  the  acoustic  wavelength  for  small  Oq,  i.e.,  1.7  mm  at  420  kHz  and  3.6  mm  at 
200  kHz.  However,  other,  corner-reflecting  features  such  as  larger  pebbles  and  rocks 
could  very  well  be  contributing  to  the  scattering  in  a  major  way,  and  it  is  unlikely 
that  Bragg  scattering  alone  is  the  source  of  riverbed  scattering. 

Finally,  radar  observations  [13]  of  backscattering  fi:om  gravel  roads  at  the  equiv¬ 
alent  radar  wavelength  in  air  also  show  the  scattering  to  increase  with  frequency  as 
we  observe. 
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3  MEASUREMENTS  IN  THE  WOOD  RIVER 

3.1  Site  Description  and  Geometry 

The  Wood  River  measurements  were  made  on  4  and  5  July  1999  at  a  site  on  the 
right  bank  about  18  miles  upstream  from  the  river’s  entry  into  Nushagak  Bay  (Fig. 
10).  This  site  represents  an  interesting  contrast  to  the  one  on  the  Kenai  River.  First, 
the  Wood  River  originates  from  an  outlet  of  Lake  Aleknagik.  The  water  is  remarkably 
clear  and  free  of  turbidity,  and  for  nearly  50  years  ADF&G  has  been  using  this  site, 
and  a  similar  one  on  the  left  bank,  to  visually  estimate  escapements  of  migrating 
sockeye  salmon  {Oncorhyncus  nerkd).  Visual  counts  of  salmon  are  made  from  atop 
20-ft  counting  towers  placed  near  the  banks  at  distances  to  the  water  exceeding  30  ft. 
Second,  the  overall  geometric  scale  in  terms  of  sonar  ranges  and  depths  is  smaller  at 
the  Wood  River  site.  Currents  in  the  Wood  River  are  considerably  faster,  typically 
1-1.5  m/s,  and  the  smaller  sockeye  salmon  tend  to  swim  in  tight  bands  1-2  m  wide 
within  1-5  m  of  the  river  bank  [14]  to  take  advantage  of  more  turbulent  flow  conditions 
and  thereby  conserve  energy.  Thus,  our  acoustic  measurements  at  this  site  were 
limited  to  ranges  <6  m  (with  a  few  measurements  made  beyond  this  range)  and  to 
depths  less  than  about  1  m. 


Figure  10.  Map  of  the  Wood  River  showing  the  location  of  the  sonar  test  site. 
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Figure  11  shows  the  position  of  the  HTI  transducer,  with  Ht  =  0.19  m,  and  the 
bottom  contour  out  to  about  5  m  in  range,  which  was  the  limit  that  could  be  measured 
directly  with  a  probe.  The  Biosonics  location  was  similar,  with  Ht  =  0.20  m.  The 
bottom  slope  varied  slowly  with  range,  with  an  average  slope,  $s,  of  8°.  No  CTD 
measurements  were  made  at  this  site  as  the  water  was  considered  both  well  mixed  and 
free  of  salinity.  (There  is  a  small  tidal  effect  that  influences  the  river  current,  which 
we  discuss  later.)  The  water  temperature  was  8°C,  setting  the  sound  speed,  c,  equal 
to  1438.5  m/s,  and  we  again  use  a  =  0.013  dB/m  at  200  kHz  and  a  =  0.060  dB/m 
at  420  kHz. 


Figure  11.  Experimental  geometry  for  the  acoustic  reverberation  measurements 
made  in  the  Wood  River  and  key  geometric  variables.  View  represents  a  portion 
of  the  river  cross  section  at  its  right  bank  extending  out  about  5  m  corresponding 
to  the  location  of  the  HTI  system.  The  transducer’s  position  is  shown  by  the  large 
triangle  at  range  =  0  m.  The  bottom  contour  has  been  constructed  from  depth 
readings  made  with  a  probe. 


3.2  Acoustic  Measurements 

For  the  420-kHz  system  the  optimum  transducer  pitch  angle  for  counting  fish, 
$Pj,,  was  found  to  be  8°,  and  measurements  were  made  at  this  angle  and  at  6°,  10°, 
0°,  and  —2°.  For  the  200-kHz  system  9pp.  was  found  to  be  7°,  and  measurements 
were  made  at  this  angle  and  at  5°,  9°,  0°,  and  -2°.  (The  transducer  pitch  angles 
reported  here  are  to  the  nearest  degree.)  As  at  the  Kenai  River  site,  the  optimum 
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transducer  pitch  angle  also  comes  close  to  the  average  slope  of  the  riverbed.  Key  sonar 
system  parameters  used  during  the  Wood  River  measurement  phase  are  summarized 
in  Table  2. 


Table  2.  Sonar  system  parameters  used  for  acoustic  measurements  in  the 
Wood  River. 


System 

Frequency 

(kHz) 

Source  Level 
(dB) 

0VK 

(deg) 

4>n 

(deg) 

101og,„i/) 

(dB) 

Biosonics  (circ.) 

420 

212.4 

4.7 

4.7 

-24.2 

Biosonics  (ellip.) 

420 

217.5 

7.1 

2.6 

-25.0 

HTI  (circ.) 

200 

216.3 

6 

6 

-21.3 

HTI  (ellip.) 

200 

213.8 

6 

10 

-20.2 

Figure  12  shows  the  results  of  two  ensemble-averaged  RL  measurements  (gray 
lines)  versus  range  made  with  the  200-kHz  sonar,  for  which  9p  was  0°.  The  upper 
curves  represent  measurements  made  on  4  July  at  approximately  1700,  and  the  lower 
curves  represent  measurements  made  on  5  July  at  approximately  1200.  The  models 
for  mean  RL  (black  lines)  result  from  evaluating  the  sonar  equation  using  the  data 
shown  in  Table  2  for  the  200-kHz  circular  transducer;  the  upper  curve  is  based  on 
Sy  =  —55  dB,  and  lower  curve  on  Sy  =  —66  dB. 


Figure  12.  Mean  RL  measured  at  200  kHz  versus  range  based  on  a  100-ping  en¬ 
semble  average  (gray  lines),  and  a  model  for  RL  (black  lines).  Upper  measurements 
were  made  on  4  July  at  approximately  1700,  and  the  lower  measurements  were  made 
on  5  July  at  approximately  1200. 
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Figure  13  shows  the  results  of  two  ensemble-averaged  RL  measurements  versus 
range  made  with  the  420-kHz  sonar  for  9p  =  0°  (gray  line)  and  Op  =  8°  (dashed 
line),  which  is  discussed  further  below.  The  measurements  were  made  on  4  July  at 
approximately  1700.  The  model  for  the  mean  RL  (thick,  black  line)  is  the  result 
of  evaluating  the  sonar  equation  using  the  data  shown  in  Table  2  for  the  420-kHz 
elliptical  transducer,  based  on  Sy  =  —42  dB. 
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Figure  13.  Mean  RL  measured  at  420  kHz  versus  range  based  on  a  100-ping 
ensemble  average  for  Op  =  0°  (thick,  gray  line)  and  Op  =  8°  (dashed  line),  and  a 
model  for  RL  (thick,  black  line). 

The  4  July  estimates  of  Sy  at  200  kHz  and  420  kHz  differ  by  about  13  dB,  which 
again  strongly  suggests  that  Rayleigh  scattering  is  in  effect.  The  source  of  this  scat¬ 
tering  in  the  clear  and  fast  moving  water  of  the  Wood  River  is  most  certainly  small 
bubbles  with  radii  <  10  //m  (see  appendix).  There  is  a  difference  of  11  dB  between 
the  200-kHz  measurements  made  on  4  July  and  5  July,  which  we  believe  is  linked  to 
the  tidal  cycle  that  determines  the  river  currents.  Unfortunately,  scheduled  current 
measurements  had  to  be  cut  back  owing  to  instrument  difficulties,  and  no  current 
measurements  were  made  close  to  the  time  of  the  acoustic  measurements.  However, 
the  4  July  acoustic  measurements  were  made  about  an  hour  after  low  tide  as  mea¬ 
sured  at  Clark’s  Point  near  the  river’s  entry  into  Nushagak  Bay,  while  the  5  July 
measurements  were  made  2-3  hours  after  a  20-ft  high  tide  at  Clark’s  Point.  The 
relation  between  the  tidal  elevation  at  Clark’s  Point  and  river  currents  is  such  that 
the  river  current  at  the  measurement  site  would  have  been  swifter  during  the  4  July 
measurements  [15].  The  swifter  current  would  provide  more  opportunity  for  bubble 
creation,  which  would  increase  the  acoustic  scattering  as  was  observed. 

Figure  14  is  a  sample  histogram  of  the  RL  at  a  range  of  5  m,  computed  for 
each  ping  of  the  100-ping  420-kHz  measurement  shown  in  Figure  13,  compared  with 
the  PDF  given  by  Eq.  (5)  (curve).  The  probability  of  the  statistic  in  this  case 
is  greater  than  0.5,  suggesting  that  the  PDF  model  is  also  a  good  one  for  the  RL 
(and  Sy)  observed  in  the  Wood  River,  provided  stationarity  assumptions  are  satisfied. 
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Figure  14.  Sample  histogram  of  the  ping-by-ping  420-kHz  reverberation  level  (in 
dB)  corresponding  to  the  data  shown  at  a  range  of  5  m  in  Figure  13.  Curve  is  a 
model  PDF,  and  the  probability  statistic  is  greater  than  0.5. 

For  both  the  200-kHz  and  420-kHz  measurements  made  on  4  July  between  1630  and 
1730,  the  variation  in  Sy  (±2  dB)  was  smaller  than  the  variation  of  ±4  dB  for  the 
Kenai  River  measurements.  But  the  large  change  in  Sy  between  the  measurements 
made  on  4  July  and  5  July  at  200  kHz  shows  the  degree  to  which  tidal-induced  changes 
in  river  currents,  and  therefore  bubble  production,  can  affect  Sy.  For  a  measure  of 
the  expected  spread  in  Sy  for  the  Wood  River  site,  we  thus  take  the  entire  range,  i.e., 
—55  dB  to  —66  dB  for  200  kHz  and  —42  dB  to  —53  dB  at  420  kHz.  Note  that  we 
infer  the  lower  limit  (—53  dB)  at  420  kHz  from  that  expected  for  Rayleigh  scattering, 
given  a  measured  Sy  of  —66  dB  at  200  kHz. 

Finally,  the  July  measurements  at  Wood  River  fell  well  within  the  peak  migra¬ 
tory  period  of  sockeye  salmon  and,  unlike  during  the  Kenai  measurements,  it  was 
difficult  to  find  time  segments  during  which  there  were  significant  gaps  in  fish  pas¬ 
sage.  Although  there  were  attempts  to  divert  fish  away  from  the  acoustic  beams, 
most  measurements  made  at  pitch  angles  corresponding  to  6pp  and  greater  showed 
evidence  of  scattering  from  fish.  Thus  we  were  unable  to  obtain  reliable  estimates 
of  scattering  from  the  riverbed  without  significant  contamination  by  scattering  from 
fish.  A  reasonable  starting  model,  however,  would  be  the  model  equations  for  Si,  at 
200  and  420  kHz  obtained  for  the  Kenai  measurements.  In  some  cases,  though,  pings 
showing  evidence  of  fish  scattering  could  be  removed  from  the  100-ping  ensemble  col¬ 
lected  by  the  Biosonics  system  before  computing  the  ensemble  average.  An  example 
of  such  a  correction  is  shown  by  the  dashed  line  in  Figure  13  representing  the  RL 
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measured  when  6p  was  8°,  the  optimum  pitch  angle  for  fish  detection.  The  reasonable 
agreement  with  the  modeled  RL  based  on  volume  scattering  alone  suggests  that,  as 
for  the  Kenai  River,  background  reverberation  at  the  Wood  River  site  is  set  by  Sy  for 
sonar  fish-counting  geometries. 
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4  SUMMARY  AND  CONCLUSIONS 


Environmental  acoustic  measurements  were  made  in  two  distinct  Alaskan  river¬ 
ine  environments  during  the  summer  of  1999:  in  the  slower-flowing  and  more  turbid 
Kenai  River,  and  in  the  faster-flowing  and  clear  Wood  River.  The  goal  of  these  mea¬ 
surements  was  to  understand,  and  model,  the  sources  and  magnitude  of  background 
acoustic  reverberation,  defined  as  the  average  acoustic  intensity  versus  time  (range) 
in  the  absence  of  any  contribution  by  scattering  from  fish.  The  measurements  were 
made  using  a  200-kHz  and  a  420-kHz  sonar  system,  in  conditions  and  geometries  that 
simulated  the  optimum  configuration  for  riverine-sonar  fish  counting.  In  addition, 
the  configurations  were  altered  by  changing  the  transducer  pitch  angle,  9p/m  order 
to  isolate  the  contributions  of  scattering  from  the  water  volume  and  scattering  from 
the  riverbed  to  the  total  background  reverberation. 

The  measurements  produced  estimates  of  the  calibrated,  ensemble-averaged  back¬ 
ground  reverberation  level,  RL  (in  dB  re  /iPa),  based  on  a  100-ping  ensemble  and 
estimates  of  the  variance  of  this  quantity.  We  modeled  RL  as  the  incoherent  sum 
of  two  components.  One,  is  volume  scattering  from  entrained  bubbles  and  silt 
particles,  for  which  the  key  determinate  variable  is  volume  scattering  strength  (in 
dB).  The  other,  RL^,,  is  scattering  from  the  riverbed,  for  which  the  key  determinate 
variable  is  the  riverbed  scattering  strength  Sb  (in  dB).  Note  that  a  third  contribution 
is  scattering  from  the  water  surface.  However,  this  contribution  is  negligible  in  com¬ 
parison  to,  and  cannot  be  distinguished  from,  volume  scattering  from  the  region  near 
the  surface. 

Our  study  shows  that  volume  scattering,  and  thus  Sy,  is  the  more  important 
of  the  two  primary  scattering  contributions.  For  example,  when  the  sonar  beams 
were  optimally  aimed  for  fish  counting — typically  a  few  degrees  below  horizontal  at  a 
pitch  angle.  Op,  close  to  the  average  slope  of  the  riverbed,  9s,  the  measured  RL  was 
well  modeled  by  an  equation  for  RLy  in  which  the  inputs  were  the  particular  sonar 
parameters  used  during  the  measurements  of  RL  and  Sy.  For  the  Kenai  River  site, 
we  estimated  Sy  to  be  —53  dB  ±4  dB  at  420  kHz  and  —67  dB  ±4  dB  at  200  kHz.  The 
difference,  40  log^o  420/200,  is  that  predicted  by  Rayleigh  scattering.  For  the  Wood 
River  site,  we  estimated  that  Sy  ranged  from  -55  dB  to  -66  dB  at  200  kHz  and 
from  —42  dB  to  —53  dB  at  420  kHz.  This  larger  range  was  associated  with  tidally 
influenced  river  currents,  which  would  modulate  the  bubble  density  by  producing 
more  bubbles  when  the  water  was  swifter. 

We  also  modeled  fluctuations  in  RL,  and  demonstrated  that  under  stationary 
conditions  fluctuations  at  a  fixed  range  are  well  described  by  a  probability  density 
function  (PDF)  given  by  Eq.  (5).  This  also  means  that  fluctuations  in  Sy  at  a  fixed 
range  are  described  by  the  same  PDF  when  stationarity  applies.  However,  as  noted 
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by  the  reported  range  in  estimates  of  Sy,  stationarity  can  typically  be  assumed  to 
apply  for  only  short  periods.  In  future  work  we  will  investigate  an  alternative  PDF 
that  would  apply  to  a  slowly  varying  Sy,  as  discussed  by  Dahl  and  Plant  [12], 

A  model  for  Sy  due  to  scattering  from  glacial  silt  particles  suspended  in  the  water, 
and  thus  applicable  to  the  Kenai  River,  is  presented  in  the  appendix.  Although 
the  model  estimates  of  Sy  agreed  quite  well  with  the  Sy  measured  at  200  kHz  and 
420  kHz,  we  caution  that  a  possible  added  contribution  due  to  scattering  from  bubbles 
is  difficult  to  distinguish  without  additional  analysis.  The  model  should  motivate, 
however,  further  study  on  the  relationship  between  Sy  and  turbidity,  as  measured  in 
nephelometric  turbidity  units  (NTU),  and  it  may  be  particularly  useful  for  acoustical- 
limnological  studies  in  Alaskan  lakes  containing  glacial  meltwater. 

Our  study  quantifies  how  fish  detection  using  horizontal-looking  sonars  in  rivers 
is  reverberation  limited  and  not  noise  hmited.  This  agrees  with  Trevorrow’s  study 
[16]  of  the  background  reverberation  in  horizontal-looking  sonars  in  Canada’s  Fraser 
River.  Increasing  the  source  level  {SL),  then,  increases  the  amplitude  of  not  only 
the  fish  echoes  but  also  the  background  RL,  and  thus  the  ratio  of  fish-to-background 
amplitude  remains  the  same. 

Knowing  Sy,  and  proceeding  with  the  idea  that  volume  scattering  alone  establishes 
the  background  reverberation  level,  then  makes  a  number  of  important  estimates  pos¬ 
sible.  One  is  the  minimum  fish  target  strength,  TSmin,  that  can  be  reliably  detected. 
Assuming,  for  example,  that  TSaan  must  be  at  least  10  dB  greater  than  the  equiva¬ 
lent  target  strength  of  the  volume  of  water  ensonfied,  TSyoi,  then  the  simple  relation 
follows  that  TS jnin  ^  TSyoi  +  10  dB,  where 


TSyol  =  10  logio  ( (13) 

Equation  (13)  can  be  used  to  evaluate  the  effect  of  different  combinations  of  pulse 
length,  r,  sonar  beams  as  parameterized  by  tp,  and  ranges,  R,  on  fish-detection  goals 
and  algorithms.  More  realism  can  be  introduced  in  Eq.  (13)  by  applying  statistical 
models  for  both  Sy  and  fish  target  strength.  In  terms  of  Sy  the  PDFs  discussed  here 
and  in  Ref.  12  will  be  useful.  Regarding  fish  target  strength,  Dahl  and  Mathisen  [17] 
have  shown  that  measurements  of  the  target  strength  of  salmon  at  side  aspect,  under 
controlled  conditions,  are  also  distributed  according  to  the  PDF  given  by  Eq.  (5), 
and  an  aspect  dependence  in  the  mean  TS  is  also  discussed.  More  recent  studies  by 
Burwen  and  Fleischman  [18]  and  Kubecka  [19]  address  how  fish  aspect  and  orientation 
in  the  riverine  environment  influence  TS  and  the  PDF  for  fish  target  strength. 

Another  important  estimate  that  we  can  make  based  on  knowledge  of  the  back¬ 
ground  reverberation  level  is  the  rms  error  in  angular  location  estimates,  ^rmsj  made 
with  split-beam  sonars  for  the  purpose  of  fish  tracking.  This  error  in  radians  is  given 
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by  [20] 


1 

WSNRd/A  ’ 


(14) 


where  d  is  the  effective  horizontal  separation  between  the  transducer  half  beams  in  a 
split-beam  system,  A  here  stands  for  acoustic  wavelength,  and  SNR  is  the  full  beam 
signal-to-noise  ratio.  For  this  case,  the  applicable  SNR  can  be  estimated  by 


2^‘5vol)/10 


(15) 


where  T5fish  is  the  nominal  target  strength  of  a  fish  at  side  aspect,  and  TS^oi  is  given 
by  Eq.  (13). 
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APPENDIX 

Comparative  Measures  of  Scattering  from  Bubbles  and  Glacial  Silt 
Particles  and  a  Model  Relating  Turbidity  Measures  to  Volume  Scattering 

from  Glacial  Silt  Particles 

Two  primary  sources  of  volume  backscattering  were  observed  in  this  study.  One 
is  small  bubbles  whose  radii  are  less  than  about  10  fim,  as  we  discuss  further  below. 
Bubble  scattering  would  be  in  effect  in  both  the  Kenai  and  Wood  rivers,  but  we 
would  expect  the  bubble  concentration  to  be  greater  in  the  faster  flowing,  and  more 
turbulent.  Wood  River.  The  other  source  is  particles  associated  with  glacial  silt  whose 
radii  are  also  about  10  ^lTa.  This  source  would  apply  only  to  the  Kenai  River,  which 
carries  a  load  of  glacial  silt,  as  the  water  at  the  Wood  River  site  is  very  clear  with 
negligible  turbidity. 

Figure  A1  shows  the  backscattering  cross  section,  ai,s  (expressed  in  dB),  for  a 
single  bubble  or  particle  of  glacial  silt  versus  bubble  or  particle  radius,  a,  computed 
using  equations  discussed  in  standard  texts,  e.g..  Ref.  21.  For  the  bubble  we  use 

=  l(/„//)^  -  1]2  + 


Figure  Al.  Backscattering  cross  section,  a^s  (expressed  in  dB),  for  a  single  bubble 
or  particle  of  glacial  silt  versus  bubble  or  particle  radius,  a,  for  acoustic  frequencies 
of  420  kHz  (solid  line)  and  200  kHz  (dashed  line). 
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where  /  is  acoustic  frequency,  5  is  the  total  damping  coefficient  (approximated  by 
S  =  2.55  X 10“^/^/^),  and  all  values  are  expressed  in  MKS  units.  Scattering  is  maximal 
at  a  frequency  equal  to  the  resonant  frequency,  /h,  for  a  bubble  of  radius  a^-  Since 
we  are  dealing  with  very  shallow  depths,  /r  and  ur  are  related  by  ur  3.25//r 
(MKS).  The  peaks  shown  in  the  figure  are  associated  with  resonant  scattering  and 
correspond  to  a  radius,  a,  of  about  7.7  /xm  for  420  kHz  and  16  /xm  for  200  kHz. 

For  a  silt  particle  we  use 


(Tbs  =  {kaY 


e-1  g-1 

3e  "^2^  +  1 


l2 


(A2) 


where  k  is  the  acoustic  wavenumber  (1829  radians/m  at  420  kHz  and  871  at  200  kHz), 
e  is  the  ratio  of  the  elasticity  of  the  particle  to  that  of  the  water,  and  g  is  the  ratio 
of  the  density  of  the  particle  to  that  of  the  water.  For  e  and  g  we  require  the  density 
and  compressional  wave  speed  of  a  silt  particle,  values  that  depend  on  the  specific 
mineral  properties  and  degree  of  homogeneity  of  silt.  A  reasonable  starting  point  is 
to  assume  that  the  silt  material  properties  are  like  those  of  sand  and  quartz,  with 
a  density  and  compressional  wave  speed  in  the  neighborhood  of  2650  kg/m^  and 
5100  m/s,  respectively  [22].  These  values  set  e  equal  to  32.78  and  g  equal  to  2.65. 

Our  data  are  very  clear  insofar  as  we  consistently  observed  a  difference  of  approx¬ 
imately  13  dB  in  the  estimates  of  Sy  derived  from  measurements  made  at  200  kHz 
and  420  kHz.  This  implies  that  Rayleigh  scattering,  in  which  the  scattered  intensity 
goes  as  the  fourth  power  of  frequency,  is  in  effect.  As  shown  by  the  upper  two  parallel 
lines  separated  by  12.9  dB  in  Figure  Al,  Rayleigh  scattering  applies  for  bubbles  whose 
radii  are  less  than  about  7  /xm.  Rayleigh  scattering  from  silt  applies  for  particle  radii 
up  to  at  least  200  fxm  (at  which  point  we  end  our  calculation).  Another  important 
consequence  of  Rayleigh  scattering  is  that  the  scattered  intensity  goes  as  the  sixth 
power  of  the  radius,  a.  The  abs  for  a  20-fxm  silt  particle,  for  example,  is  18  dB  greater 
than  that  for  a  10-/im  silt  particle. 

Recall  that  our  observable,  however,  is  not  abs  but  rather  the  backscattering  cross 
section  per  unit  volume,  Sy  (or  its  decibel  equivalent  Sy),  which  is  defined  by  the 
integral 

Sy=  abs{(x)N{a)da.  {A3) 

Ja 

The  distribution  function  N{a)  gives  the  number  of  particles,  or  bubbles,  per  cubic 
meter  with  radii  between  a  and  a  +  da,  and  thus  is  analogous  to  a  probability  density 
function  (PDF)  for  bubble  or  particle  radius. 

Bubbles  within  a  wide  size  range  will,  in  fact,  be  produced  by  the  breaking  of 
the  water  surface  existing  in  areas  with  fast  currents  or  areas  with  flow  obstructions 
caused  by  large  rocks.  Larger  bubbles  quickly  rise  to  the  surface,  but  smaller  ones  are 


A2  TR  2001 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


easily  advected  a  few  meters  below  the  surface  by  turbulent  diffusion.  The  persistence 
time  of  smaller  bubbles  (radii  <  10  yum)  is  lengthened  in  the  cold,  highly-oxygenated 
riverine  water,  which  invariably  contains  surfactants  associated  with  biological  activ¬ 
ity  [23]. 

Thus  we  postulate  that  bubbles  with  radii  <  10  fim  are  the  only  source  of  the 
volume  scattering  observed  in  the  shallower,  faster  flowing,  and  more  turbulent  Wood 
River.  We  do  not  have  sufficient  information  to  characterize  the  size  distribution  of 
these  bubbles  further,  other  than  the  one  linked  to  Rayleigh  scattering.  Bubbles 
are  also  likely  to  be  a  source  of  volume  scattering  in  the  Kenai  River.  We  expect 
bubble  concentration  in  the  Kenai  River  to  be  less  than  that  found  in  the  Wood 
River,  however,  since  the  current  is  slower  in  the  Kenai.  The  exception  is  wakes  from 
boat  traffic  in  the  Kenai,  which  we  have  shown  to  produce  an  enormous  increase  in 
Sy  associated  with  bubble  scattering. 

The  significant  load  of  glacial  silt  carried  in  the  Kenai  River  and  the  resulting 
turbidity,  however,  mean  that  scattering  from  suspended  particles  also  contributes  to 
the  volume  backscatter.  It  is  not  possible  to  determine  the  exact  relative  contribution 
of  each  source,  bubbles  and  silt  particulates,  to  acoustic  volume  scattering  because 
we  do  not  have  the  requisite  ancillary  data  on  in  situ  bubble  and  silt  particulate  size 
distribution  and  concentration.  We  can,  however,  make  a  first-cut  evaluation.  For 
this,  we  shall  use  some  empirical  relations  for  three  limnology  variables  that  have 
been  measured  by  scientists  from  ADF&G  in  Alaskan  lake  and  riverine  systems  that 
are  influenced  by  glacial  meltwater. 

The  first  variable  is  the  depth,  SD  (in  meters),  to  which  a  secchi  disk  can  be  seen 
(the  secchi  transparency  depth),  which  is  a  measure  of  water  transparency.  ADF&G 
biologists  recorded  SD  about  3  times  daily  on  the  Kenai  River  at  a  site  near  our  exper¬ 
imental  area.  From  these  recordings  [24]  we  estimate  that  the  SD  was  approximately 
0.6  m  at  the  time  of  our  volume  scattering  measurements. 

The  second  variable  is  turbidity  expressed  in  nephelometric  turbidity  units,  or 
NTU.  We  will  see  that  NTU  is  a  key  variable  for  examining  acoustic  volume  scattering 
from  glacial  silt,  because  of  its  strong  correlation  to  both  particle  concentration  and 
mean  particle  radius  [25].  Turbidity  levels  are  inversely  proportional  to  the  light  com¬ 
pensation  depth  as  determined  by  a  secchi  disk  reading  [25],  and  an  empirical  relation 
predicting  SD  as  a,  function  of  turbidity  expressed  in  NTU  has  been  demonstrated  by 
Koenings  and  Edmundson  [26]  using  data  from  several  Alaskan  lakes  influenced  by 
glacier  meltwater.  Edmundson  [27]  kindly  re-examined  data  from  15  Alaskan  lakes 
with  turbidities  that  encompass  the  range  of  measurements  in  the  Kenai  River  and 
derived  the  following  relation  enabling  us  to  use  SD  as  a  model  predictor  of  NTU: 


logio(NTU)  =  0.952  -  1.093  logio(SD) 


(A4) 
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with  squared  correlation  coefficient  =  0.89.  Using  Eq.  (A4)  and  SD  =  0.6  m  gives 
an  estimate  of  NTU  =  16. 

The  third  variable  is  mean  particle  size,  PS  (or  diameter).  Edmundson  and 
Koenings  [25]  also  show  that  PS  is  also  strongly  correlated  with  NTU.  We  have  re¬ 
digitized  the  data  presented  in  Figure  lA  of  Ref.  25  in  order  use  NTU  as  a  model 
predictor  of  PS  in  microns.  The  result  is 

log, q{PS)  =  -0.355 \og,o{NTU)  +  1.704  ,  (A5) 

with  =  0.75.  Using  Eq.  (A5)  with  NTU  =  16  puts  PS  equal  to  about  19  /xm,  for 
a  mean  particle  radius,  a,  between  9  and  10  ^m. 

To  proceed  further,  we  require  a  description  of  N{a).  The  primary  author  (Dahl) 
has  examined  an  unpublished  set  of  histograms  of  particle  sizes  in  samples  of  glacial 
meltwater  taken  in  the  mid-1980s  by  scientists  from  ADF&G.  A  Rayleigh  probability 
density  function  (PDF)  appears  to  be  a  reasonable  model  for  these  histograms.  The 
Rayleigh  PDF  as  a  function  of  particle  radius,  a,  is 

Paia)  =  (A6) 

The  PDF  is  governed  by  aa,  which  is  related  to  the  expected  value  of  particle  radius, 
E{a),  by  E{a)  =  By  hypothesis,  we  equate  E{a)  to  mean  particle  radius, 

or  PS/2,  as  determined  by  Eq.  (A5).  The  particle  size  distribution,  N{a),  used  in 
Eq.  (A3)  is  PN  x  Pa{o.),  where  PN  is  the  total  number  of  particles  per  cubic  meter. 
We  can  now  estimate  PN  using  Eq.  (A3),  along  with  our  nominal  observed  values  of 
Sy  for  the  Kenai  River  which  are  —67  dB  at  200  kHz  and  —53  dB  at  420  kHz.  These 
values  for  Sy,  combined  with  the  above  assumptions,  give  a  PN  on  the  order  of  10^^ 
silt  particles  per  cubic  meter. 

This  is  plausible  value  for  PN.  However,  there  are  still  many  issues  to  resolve. 
One  is  the  degree  of  flocculence  of  the  suspended  silt  particulates,  i.e.,  the  degree  to 
which  the  particulates  have  coalesced.  Another  is  the  relative  contribution  of  bubbles 
to  the  total  volume  scattering. 

Thus,  the  main  purpose  of  this  analysis  is  to  serve  as  a  guide  and  to  motivate 
further  study  of  the  relationship  between  Sy  and  turbidity  as  measured  in  units  of 
NTU.  For  example,  a  longer  time  series  of  simultaneous  measurement  of  Sy,  NTU,  and 
SD  at  the  minimum,  and  if  resources  permit,  histograms  of  particle  size  distribution, 
would  help  clarify  several  issues. 
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